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LINEARIZATION COEFFICIENTS FOR ORTHOGONAL 
POLYNOMIALS USING STOCHASTIC PROCESSES 1 

By Michael Anshelevich 

University of California, Riverside 

Given a basis for a polynomial ring, the coefficients in the expan- 
sion of a product of some of its elements in terms of this basis are 
called linearization coefficients. These coefficients have combinatorial 
significance for many classical families of orthogonal polynomials. 
Starting with a stochastic process and using the stochastic measures 
machinery introduced by Rota and Wallstrom, we calculate and give 
an interpretation of linearization coefficients for a number of polyno- 
mial families. The processes involved may have independent, freely in- 
dependent or q-independent increments. The use of noncommutative 
stochastic processes extends the range of applications significantly, 
allowing us to treat Hermite, Charlier, Chebyshev, free Charlier and 
Rogers and continuous big q-Hermite polynomials. 

We also show that the q-Poisson process is a Markov process. 

1. Introduction. Let {P n } be a family of polynomials orthogonal with 
respect to a measure /i on the real line. One standard combinatorial question 
is to calculate the moments of the measure, m n = (x n ), where we denote by 
(•) the integral (expectation) with respect to \x. For many classical families 
of polynomials these moments are positive integers or, more generally, poly- 
nomials in parameters with positive integer coefficients. These coefficients 
beg a combinatorial interpretation, and there exists a large body of work to 
this effect. 

A more general question one can ask is to calculate the linearization co- 
efficients. That is, for (ni,n2, . . . , we are interested in the expectations 
(P ni Pn 2 ■ ■ ■ Pn k ) • The name stems from the fact that these are the coefficients 
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in the expansion of products of this type in the basis {P n }, that is, expan- 
sions as sums of orthogonal polynomials. Again, many of these coefficients 
are positive integers, and so they "count something." 

A combinatorial approach to this problem is to construct explicit bijec- 
tions between structures counted by the linearization coefficients and struc- 
tures of known cardinality; see, for example, [10]. In this paper we take a 
different route, and consider a probabilistic interpretation of certain coef- 
ficients. The connection to combinatorics is provided by the fact that the 
moments of a measure are sums, over all set partitions, of products of cu- 
mulants of that measure. We will see that certain linearization coefficients 
can by described in a similar way. The machinery we use is that of stochas- 
tic measures, first introduced by Rota and Wallstrom in [12]. In a number 
of previous papers we extended this machinery from the usual to the non- 
commutative stochastic processes. This extends the number of polynomial 
families that we can handle, and so we not only obtain a nice interpretation 
of known results, but some new results as well. In particular, we show that 
the linearization coefficients for the continuous big q-Hermite polynomials 
([11], 3.18) are based on the number of the inhomogeneous set partitions, 
with an extra statistic counting the number of "restricted crossings" of such 
partitions. 

To be more specific, for each family of polynomials in this paper and 
the related family of measures of orthogonality, we introduce a, possibly 
noncommutative, stochastic process {X(t)}. Then for this process, we in- 
troduce a further family {4>k(t)} of other stochastic processes, which we call 
full stochastic measures. These objects are orthogonal, and have clean lin- 
earization formulas. On the rare occasions when these objects are polynomi- 
als in the original process X, these formulas translate into the linearization 
formulas for polynomials. 

Another property, which always holds for the full stochastic measures 
and which in these cases is shared by the orthogonal polynomials, is the 
martingale property. The Markov property for the g-Brownian motion was 
shown in [6] using the Gaussian properties of the process. Using the above 
fact for the stochastic measures, we show that the g-Poisson process is also 
a Markov process. 

The paper is organized as follows. In Section 2 we describe general combi- 
natorial properties of combinatorial stochastic measures. Section 3 is based 
on the results of [12] about processes with independent increments, and gives 
the linearization coefficients for the Hermite and Charlier polynomials. Sec- 
tion 4 is based on the results of [1, 3] about processes with freely independent 
increments, and gives the linearization coefficients for the Chebyshev poly- 
nomials of the 2nd kind and the free Charlier polynomials. Section 5 is based 
on the results of [2] about g-Levy processes, and gives the linearization co- 
efficients for the continuous and continuous big g-Hermite polynomials. It 



LINEARIZATION COEFFICIENTS 



3 



also requires some new results about the g-Poisson process. The proofs of 
these results are contained in the Appendix, which also contains the proof 
of the Markov property for the g-Poisson process. 

2. Combinatorial stochastic measures. Let (A, E[-]) be a noncommuta- 
tive probability space. That is, A is a finite von Neumann algebra, and E[-] 
is a faithful normal tracial state on it. The commutative case is included 
in this setting when A = L°°(A,P) for A a measure space, P a probabil- 
ity measure, and E[-] the expectation with respect to P. Let {X(t)} be an 
operator-valued stochastic process whose increments are stationary with re- 
spect to the state E[-] and independent in a certain sense; see Sections 3-5 
for examples of such conditions. Denote by V{n) the collection of all set 
partitions of a set of n elements. For a set partition ir = {Bi,B>2, ■ ■ ■ ,Bi}, 
temporarily denote by c(i) the index of the class B c ^ to which i belongs. 
Then the stochastic measure corresponding to the partition tt is 




dX(s cW )dX(s c{2 )) ■ • • dX(s cijl )). 



all Si's distinct 

In particular, denote by A n = Stj the higher diagonal measures of the pro- 
cess defined by 

A n (t) = ( (dX(s)) n , 
and by ip n = Stg, the full stochastic measures defined by 

ipn(t)=[ [M) n dX(si)dX(s 2 ) ■ ■ -dX(s n ). 

all Si's distinct 

Here the integrals are defined by approximation, as follows. 

Definition 2.1. Let 1 = {Ii}fLi be a subdivision of the interval [0,t) 
into disjoint half-open intervals h = [ai,ai+i), = a\ < a2 < ■ ■ ■ < ajv < 
a^v+i = t. Denote by 5(2) = maxj Let tt £ V(n) and {A(s)} be a (possi- 
bly noncommutative) stochastic process. Define 

TO = {it 6 {1, 2, . . . , N} n : u(i) = u(j) ^i^j} 

and 

n 

SW(t;J)= Y[( x ( a u(i)+i) -A(a u(i) )). 

3e[jv]»t=i 

Finally, define 

St,r(t) = lim St n (t;I) 

8{2)->0 

if the limit exists. 
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The existence of the limits has been proven under various conditions, see 
Sections 3-5 for the more precise description. For the purposes of this section 
we will assume that the limits exist and consider purely combinatorial facts. 
The most pertinent of these corresponds to linearization or, in the context 
of stochastic integration, to the Ito formula. Set n = J2j=i n j- Denote by 

the partition whose classes are intervals of consecutive integers of lengths 
n\ , ri2, ■ ■ ■ , rii-. Denote 

7 7 (ni,n 2 , • • • ,n k ) = {-rr £ V(n) :tt A 7r nij „ 2) ... jns , = 6}, 

the partitions inhomogeneous with respect to 7r ni na ... nfc , that is, the col- 
lection of all partitions which do not put together elements of the k distin- 
guished subsets in the same class. For example, 

P(2,2) = {{(1)(2)(3)(4)},{(1,3)(2)(4)},{(1,4)(2)(3)}, 

{(1)(2,3)(4)},{(1)(2,4)(3)},{(1,3)(2,4)},{(1,4)(2,3)}}. 

Then 

k 

(i) II^Mt)= E sw(*). 

3=1 ■K€V(n\,n2,...,n k ) 

For a fixed subdivision X, the statement 

k 

3=1 n£V(ni,ri2,...,n k ) 

is purely combinatorial; see [12], Theorem 4, or [1], Proposition 4, for its 
proof. The fact that the relation (1) also holds in the limit will again be 
treated in each of the subsequent sections separately. 

Denote R n (t) =E[St 7r (t)] and R n (t) = E[A n (t)]. Here R n is the nth gen- 
eralized cumulant of the process; for a process with independent increments 
it is the usual cumulant. Then 
r k 



(2) E 



I j=l 



R*(t). 

nGVfoi ,ri2,..., n-k) 



For a centered process, R\ = 0. In all examples we will consider, this will 
imply that R v = if ir contains a singleton class (a class consisting of one 
element). One consequence of this fact is that 

(3) E[Mt)Mt)]=o 

for n^h. That is, full stochastic measures of different orders are orthogo- 
nal. Thus, in general, we may consider the stochastic measures as analogs 
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of orthogonal polynomials, and in this case formula (2) describes their lin- 
earization coefficients. The purpose of this paper is to describe examples 
when stochastic measures are, in fact, polynomials in the original process. 
If ip n (t) = P n (X(t)), equation (3) says that the polynomials {P n } are or- 
thogonal with respect to the distribution fit of X(t) (which is a probability 
measure on M). So their linearization coefficients are precisely 



p, 



Ilk I 



E 



E 



n<M*) 

Lj=l 



Lj=l 

Moreover, in all examples below X(t) has infinite spectrum (takes on in- 
finitely many values). So if St^t) is also a polynomial in X(t), equation (1) 
implies the equality of the corresponding polynomials. 

Another property which holds for some orthogonal polynomials, but which 
always holds for stochastic measures, is the martingale property. 



Proposition 2.2. For t>0, let At be the von Neumann algebra gen- 
erated by the set {X(s) :s < t}. Assume the following: 

(a) There exist consistent conditional expectations {Et[-]} from A onto 
each At preserving the expectation E[-]. 

(b) The process {X(t}} is centered, that is, E[X(t)] = for all t. 

(c) The increments of the process are singleton independent. That is, 
given a collection of intervals [sj,tj) C R+, j = 1, 2, . . . , k such that for some 
i, 

[si,ti)n ( UMi)) =0, 

then E[(X(ti) - X( Sl )) . . . (X(U) - X( Sl )) . . . (X(t k ) - X(s k ))) = 0. 

(d) The limit defining ijj n (t;X) exists in the I? -norm with respect to E[-]. 
Then the process ip n (t;X) is a martingale with respect to the filtration 

{At}- That is, for all s <t, 

E s m;X)]=^{s;X). 



See the Appendix for the proof. 



3. Processes with independent increments. Let {X(t)} be a process 
with stationary independent increments, and, thus, a Levy process. Then 
by the results of [12], the integrals defining stochastic measures exist as 
limits in probability. Moreover, it is not hard to show that in this case for 
n = {B 1 ,B 2 ,...,B l }, 

(4) St w (t) = ip(t- A )Bl | , A| fla |, . . . , A| Bj |). 
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Here we are using a slightly more general definition of a stochastic mea- 
sure where different factors in its defining integral may come from different 
processes: 

i,(t;(X^\X^,...,X^)) = J m „ dX ( - l \s 1 )dX (2 \s 2 )---dX^(s k ). 

all Si's distinct 

See [3] for details. Throughout the paper we will consider stochastic pro- 
cesses for which the diagonal measures are affine functions in the original 
process X. Two types of processes that have this property are generalized 
Brownian motions and generalized Poisson processes. 

A stochastic measure is multiplicative if E[St 7r (t)] = Ilse^ E[Ai s i (t)] . Both 
stochastic measures in this section are multiplicative. For a multiplicative 
measure 

(5) R K (t)=*[[R ]B] (t), 
and the sum on the right-hand side of (2) is equal to 

e n 

In this case, (1) follows from [12], Theorem 4. 

Notation 3.1. Denote by Vi^in) partitions whose classes consist only 
of one or two elements, otherwise known as "matchings," and by T > 2{n), the 
collection of all pair partitions, otherwise known as "perfect matchings." 
Denote by s(tt) the number of singleton (one-element) classes of it, and by 
S2(tt), the number of two-element classes. 

3.1. Hermite. If {X(t)} is the Brownian motion, then by the strong law 
of large numbers ^(t) = t, A m (£) = for m > 2. Moreover, it follows from 
the Kailath-Segall formula (see [12], Theorem 2) that ip m (t) = H m (X(t),t). 
Here H m {x,t) are the Hermite polynomials, defined by the recursion rela- 
tions 

xH m (x, t) = H m+ i (x, t) + mtH m -\(x, t). 
It follows from (4) that for n S Pi,2, 

St n (t)=t s ^H s{w) (X(t),t) 
and they are otherwise. Therefore, (1) gives 

k 

]jH nj (x,t)= J2 t s ^H s{n) (x,t). 

3=1 TreVi t 2(ni,n2,...,n k ) 
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In particular, 

^ J] H n . (x,t)j= e' 2 \V 2 {n 1 ,n 2 ,...,n k )\. 
This formula is well known and surely quite old. Since (H m (x,t) 2 ) = m\t m , 

k n 1 

\H n .{x,t)=Y J —t in - m)/2 \V2(ni,n 2 ,...,n k ,m)\H m (x,t). 

j=l m=U 

3.2. Centered Charlier. If {X(t)} is the centered Poisson process, then 
from [12], Proposition 7, A m (i) = X{t) +t for m > 2. Moreover, ift m (t) = 
C m (X(t),t). Here C m (x,t) are the centered Charlier polynomials, defined 
by the recursion relations 

xC m (x,t) = C m+ x(x,t) +mC m (x,t) + tmC m -i(x,t). 

It follows from (4) that 



i=0 

Therefore, 

k \ir\-s(ir) 



IIc^(M)= E E ( W 7 W Wh(M)- 

i=i 7rG-P(ni,n 2 ,...,n fc ) i=0 V ' 



In particular, 



A- 



(n c ^M))= e * w - 

\j=l / ■w£V(n 1 ,n 2 ,...,n k ) 

s(tt)=0 

This formula appears in [15] and a number of later sources. Since (C m (x, t) e 
m\t m , 



n^(^*)=E^ e t^- m c m ( X ,t). 

j=l m=0 ' n£V(ni, ri2,..., n k ,m) 

s(tt)=0 



Note that the noncentered polynomials, here and in the subsequent sections, 
will have exactly the same linearization coefficients. 
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Outer classes 




i 



2 



3 




A 5 /, 



7 



9 



It) 



Inner classes 



Fig. 1. A noncrossing partition of 10 elements with 3 inner and 3 outer classes. 

4. Processes with freely independent increments. The notion of free in- 
dependence was introduced by Voiculescu [14] in the context of operator 
algebras. X,Y £ A are freely independent if, whenever 



This property is easily seen to be incompatible with, but is parallel to, the 
usual independence. Free probability is by now quite a rich theory which is 
based on this notion; see [14] for an overview. In particular, there is a well- 
developed theory of free cumulants, free infinitely divisible distributions and 
limit theorems, and processes with freely independent increments. 

In this section, let {X(t)} be a process with stationary freely independent 
increments, and, thus, a free Levy process. It was shown in [1, 3] that in this 
case the integrals defining stochastic measures exist as limits in the operator 
norm, as long as the operators {X(t)} are bounded. Moreover, it was shown 
in [1] that St,,- = unless ir is a noncrossing partition. Here a partition ir is 

7T 

noncrossing if there are no i < j < k < I with i ~ k, j ~ I, i ^ j. 

In this case, (1) follows from [1], Proposition 4. For the analog of the 
formula (4), we need a new notion. For a noncrossing partition, we distin- 
guish the classes that are inner, or covered by other classes, and outer. See 
Figure 1 for an example. 

Proposition 4.1. Let ir be a noncrossing partition with outer classes 
£?l, . . . , -B^t,-) and inner classes C\, . . . , Cu^ . Then 




i{ix) 



sww = n • ^( a ibii(*)> a ib 2 i(*)> • • • > A| B 



i=l 



LINEARIZATION COEFFICIENTS 9 
Proof. This is a particular case of the main theorem of [3]. □ 

Notation 4.2. Denote by NC{n) the lattice of noncrossing partitions, 
and by NCi )2 (n), NC 2 (n), NC(ri\, . . . ,rik), and so on, the corresponding 
subsets of NC(n). Denote by si(ir) the number of inner singletons of tc. 
Denote by o(ir) and i(ir) the number of outer and, respectively, inner classes 

of 7T. 

Free stochastic measures are not multiplicative in general. However, (5) 
does hold for n € NC{n). So for a free stochastic measure, the general lin- 
earization coefficients are 

" k 

e n<M*) = e m*)= e n%i(*)- 

-3=1 J 7reiVC(ni,ra2,...,n fc ) ireNC{n\,n2,...,n k ) B£n 

4.1. Chebyshev. There is a free version of the central limit theorem, with 
independent variables replaced by freely independent ones. The limit distri- 
bution in this theorem is the semicircular distribution. A process (consisting 
of noncommuting operators) {X(t)} with stationary freely independent in- 
crements all of which have (scaled) semicircular distributions is the free 
Brownian motion. 

If {X(t)} is the free Brownian motion, then from [1], A 2 (t) = t, A m (i) = 
for m > 2. Moreover, by [1], Corollary 8, ip m (t) = U m (X(t),t). Here U m (x,t) 
are the Chebyshev polynomials of the second kind, defined by the recursion 
relations 

xU m {x,t) = U m+ l(x,t) +tU m -l(x,t). 

It follows from Proposition 4.1 that for ir 6 NCi^(n) and si^w) = 0, 

St n (t)=t s ^U s(7T) (X(t),t) 
and they are otherwise. Therefore, by (1), 

k 

Y[U n] (x,t)= J2 t S ^U sM (x,t). 

3=1 TTeNC lt2 (ni,n 2 ,...,n k ) 
si(ir)=0 

In particular, 

^ II U n 4 {x,t)^ = t n/2 \NC 2 (n 1 ,n 2 ,...,n k )\. 

This formula has essentially appeared in [8], in a slightly different guise (they 
count the number of Dyck paths). Since (U m (x,t) 2 ) = t m , 

k n 

J] U n . (x, t)=Y, t^- m ^ 2 \NC 2 ( ni ,n 2 , ...,n k , m)\U m (x, t). 

j=l m=0 
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4.2. Centered free Charlier. The distribution of the sum on n freely in- 
dependent Bernoulli ((1 — ^)<5o + ^6\) variables converges, as n — > oo, to a 
distribution which is naturally called the free Poisson distribution with pa- 
rameter t. It is also known as the Marchenko-Pastur (or, for t=l, Wishart) 
distribution. A process with stationary freely independent increments such 
that the increments have free Poisson distributions is the free Poisson pro- 
cess. 

If {X(t)} is the centered free Poisson process, then by [1], Corollary 4, 
A m (t) = X(t) + t for m > 2. By [1], Corollary 10, ^ m {t) = C 0| m (*(*),*)■ 
Here Co ;m (x,t) are the centered free Charlier polynomials, defined by the 
recursion relations 

xC 0i o(x,t) = C 0> i(x,t), 

xC 0>m (x,t) = C , m+ i(x,t) + C 0jm (x,t) + tC 0>m -i(x,t) 

for m > 0. They are, of course, orthogonal with respect to the free Poisson 
distribution. 

It follows from Proposition 4.1 that for ir 6 NC(n) and si(n) = 0, 

o(ir)-s(n) + si(ir) , . 

st. = f« £ ( 0( ^ " s( f + s * (7r) ) t^ yi (X(t),t) 

and they are otherwise. Therefore, by (1), 

k 

]JC , nj (x,t) 



o(n)— s(ir)+si(ir) 

E E 

n€NC(m,n2,— ,n k ), 1=0 



o(ir) — s{n) + si (it) 
I 



xt lW+i C M7r) _zOM). 



In particular, 



k 



l[Co, nj (x,t))= Yl tH - 

W =1 ' ireNC(ni,n 2 ,...,n k ) 

s(n)=0 



Since (C , m (x, t) 2 } 



fc n 

I] Co^Cx.t) = £ E t^- m C , m (x,t). 

j=l m=0 Tr£NC(ni,n2,...,n k ,m) 

s(n)=0 
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5. Processes on a g-deformed full Fock space. 

5.1. q-Fock space. Consider the Hilbert space L 2 (R + ,dx). Let 

oo oo 

F alg (L 2 (R + )) = @L 2 (R + ,dx)® k = ($L 2 (R k + ,dx® k ) 

fc=0 k=0 

be its algebraic Fock space. Here the Oth component is spanned by the 
vacuum vector $7. Then (•, -)o defined by 

(/l ® ■ • ■ ® fk, 91 ® ■ • • ® 9n)o = 5kn(fl,9l) ■ ■ ■ (fk,9k) 

is an inner product on the algebraic Fock space, where (-, ■) is the standard 
inner product on I? (R + , dx) . Define the operator P q by 

P q {fl®---®fn)= E ^ (<T) /a(l)®---®/a(n), 
crgSym(n) 

where Sym(n) is the permutation group and i(a) is the number of inver- 
sions of a. According to [7], this operator is strictly positive for — 1 < q < 1. 
Denote (•, -) q = (-,P g -)o. Then this is also an inner product, and we denote 
by ^"q(L 2 (R+)) the completion of jF a i g (L 2 (]R+)) with respect to the corre- 
sponding norm, and call it the g-deformed full Fock space. 

For / G L 2 (R + ) nL°°(R + ), define creation, annihilation and preservation 
operators on the g-Fock space JF (? (L 2 (R + )) by 

a*(/)(0) = /, 

a*(f)(9i®---®9n) = /®ffl® ■■■<S>9n, 

o(/)(n) = o, 
<f)(g) = {f,g)n, 

n 

a(f)(gi ® • ■ • ® 5n) = E ^' 5ffc ) fi ' 1 ® ■ ' ■ ® 9fc ® • • • ® 5n, 
fc=i 

p(/)(n) = o, 

n 

® • • • ® 5n) = E ® 51 ® • • • ® 5fc ® • • • ® 5n, 

fe=l 

where ^ means "omit feth term." For / real-valued, is self-adjoint, and 
a(f) and a*(/) are adjoints of each other. 
The noncommutative stochastic process 

X(t)=a*(l m ) + a{l m ) 

is, by definition, the g-Brownian motion, and the process 

X(t)=a*(l m ) + a(l m )+p(l m ) 
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is the centered g-Poisson process. Let A be the von Neumann algebra gener- 
ated by {X (t)} te \ 0>oo ) > an d let E[-] be the vacuum vector state E[-] = (Q, -Q,). 
Then for q = 0, these processes are the free Brownian motion and the cen- 
tered free Poisson process, while for the degenerate case q = l, they give the 
corresponding classical processes. 

Notation 5.1. Let ir be a partition. Denote by rc(7r) the number of 
restricted crossings of ir. Here a restricted crossing of ir is a 4-tuple i < j < 
k < I such that i ~ k, j ~ I, and k = min{r :r > i,r ~ i}, I = min{r :r > 
3i r ~ j}- See Figure 2 for an example. Also, define the singleton depth sd(-7r) 
of 7r to be the sum of depths, d(i) = \{j\3a, b € Bj : a < i < b}\, over all the 
singletons (i) of it. 

For — 1 < q < 1, denote 

n-l -j _ n n 

[o] 9 = o, [n^E^-r^ and W«' = IIbV 

The stochastic measures for the g-Levy processes are described in a forth- 
coming paper [4]. However, the functionals R n are known to be well defined, 
and the following analog of (5) holds. 

Proposition 5.2 ([2], Theorem 3.8). 

R*(t) = q IcM U R \B\( t ')- 

BGtt 

Therefore, if (1) holds, the linearization coefficients are 
E 



J2 Mt)= E q rc[n) UR\B\(t). 



Proposition 5.3. If X(t) is the q-Brownian motion, then the limit 
defining Ak(t;X) exists in the L 2 -norm with respect to E[-], and equals 
A2(i) = t, Afc(t) = for k > 2. Similarly, if X(t) is the centered q-Poisson 
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process, the limit defining Ak(t;X) exists in the L? -norm with respect to 
E[), and equals A k (t) =X(t) +t for k>2. 

See the Appendix for the proof. 

5.2. q-Hermite. If {X(t)} is the q-Brownian motion, then ip m (t) = Hq tTn (X(t),t). 
Here the polynomials H qjTn (x, t) are a scaled version of the continuous (Rogers) 
g-Hermite polynomials, defined by the recursion relations 

m+ l(x,t) + t[m] q Hq 

,m—l (x,t). 

The measure of orthogonality of these polynomials, and so the distribution 
of X(t), is the most common version of the g-Gaussian measure; see [6] or 
(with a slightly different normalization) [11], 3.26. 

Lemma 5.4 ([2], Proposition 6.12). For tt <G V\,2, 

St n = q^+ sd ^t s ^Hq A7T) (X(t),t) 

and they are otherwise. Here the limit in the definition of the stochastic 
measure is in L°°~ , that is, in LP for any 1 < p < oo, with respect to E[-]. 

Therefore, 

fl H q , nj (x,t) = £ q ^ +sd ^t s ^H qAn) (x,t). 

3=1 ir€Pi,2(ni,n2,...,nk) 
In particular, 

\j=l ' Tv&P2(ni,n,2, ...,nfc) 

This formula has appeared in [9]. Since (H q ^ m (x,t) 2 ) = [m\q\t m , 

k n -I 

II H q , nj (x,t) = £ T^Tl tin - m)/2 E ^H q , m (x,t). 

3=1 "i=0 L J<? - TreV 2 {n 1 ,n2,...,n k ,m) 

5.3. Centered big q-Hermite. Let {X(t)} be the centered q-Poisson pro- 
cess. 

Proposition 5.5. For the centered q-Poisson process, if) m (t) = C q ^ m {X(t),t). 
Here C q ^ m are a scaled version of the centered continuous big q-Hermite 
polynomials, which in our context are q-analogs of the Charlier polynomials. 
They are defined by the recursion relations 

xC q , m (x,t) = Cq tm+1 (x,t) + [m] q C q , m (x ,t) + t[m]qC q , m -l(x ,t) . 

In particular, the stochastic measures ip m are well defined, with the limits 
taken in the I? -norm with respect to E[-]. 
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See the Appendix for the proof. Also, see [13] for a detailed description 
of the measure of orthogonality of these polynomials. 
From Proposition 5.2, for s(ir) = 0, 

(6) R K (t) = q r <«hW, 

and they are otherwise. 

Proof of (1). We start with the known combinatorial formula 

k 

n^ nj (t;Z)= £ St n (t;l). 

j=l ■K£'P{n 1 ,n2,...,n k ) 

By the results of [4], St 7r (t; J) converges to St 7r (t) in L 2 . By the previous 
proposition, i/j nj (t;2) converges to tpnjit) m L 2 . Thus, ip ni (t;2)tj} n2 (t;2) 
converges to t/j ni (t)ip n2 (t) in L 1 . On the other hand, it also converges to 
X] 7re p( niin2 ) St 7r (t) in L 2 . Therefore, the two latter expressions are equal, 
and ij) ni (t;2)ijj n2 (t;2) converges, in fact, in L 2 . By induction, we conclude 
that n}=i (t;T) converges in L 2 to E ff6 7»(ni,na,...,n h ) St «-(*)- □ 

We conclude that 

(n<w*,t))= e ? rc(7r) * kl - 

\j"=l / 7re"P(ni,n 2 ,...,n fe ) 

s(vr)=0 

Since (C wn (x,t) 2 > = [m] ff !t m , 

I1 C «(^) = E m E g rcW i khm O), m (*,t). 

J=l m=0 ^ J' - TTeV(ni,ri2,...,n k ,m) 

s(vr)=0 

5.4. Limiting relations. The results of the previous sections can be ob- 
tained as the limits of the results of this one. For the continuous (Rogers) 
g-Hermite polynomials, taking q = 1 gives the formulas for the Hermite poly- 
nomials, while taking q = gives the formulas for the Chebyshev polynomi- 
als. For the continuous big g-Hermite polynomials, taking q = 1 gives the 
formulas for the Charlier polynomials, while taking q = gives the formu- 
las for the free Charlier polynomials. Note that in the latter case we only 
recover the linearization coefficients themselves, not the expressions for the 
products of polynomials as sums over partitions. 

Finally, consider the process 

X(t,a) =a*(l[ ,t)) +a(l[o,t)) + ap(l[o,t))- 
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For this process, A m (t) = a m ~ 1 X(t) + a m ~ 2 t and R m (t) = a m ~ 2 t for m > 2. 
Therefore, for s(ir) = 0, 

R n = q rc ^k^a n ~ 2 ^, 

and they are otherwise. Also, for this process, ip m (t) = P gjmi0 ,(X (t) , t), 
where 

%Pq,m,a t) — Pq,m+l,a t) ~\~ Ot\jTl\qPq rna (x, t) -|- t\lTl\qPq m -.\ a (x, t) . 

We conclude that 

\j=l I ireV(ni,n 2 ,—,n k ) 

s(tt)=0 

For a = 1 , this gives the g-Poisson process and the continuous big q-Hermite 
polynomials. On the other hand, for a = 0, this gives the q-Brownian motion 
and the continuous (Rogers) g-Hermite polynomials. In the linearization 
formula, the only partitions with a nonzero contribution are those with n = 
2\ir\ and without singletons, that is, pair partitions. 



APPENDIX 

g-Levy processes. We briefly review the definition of more general q- 
Levy processes and their stochastic measures; see [2] for more details. Let 
V be a Hilbert space, and consider H = L 2 (R + ,dx) <g> V. Define F a i s (H), 
T q (H), E[-] and, for £ G H, a(£) and a*(£) as in the beginning of Section 5 
for V = C. For T an essentially self-adjoint operator on H, define p(T) on 
? q {H) by 

p(T)(fi) = 0, 

n 

P(T)($1 ® • • • ® &») = E ff*" 1 ^*) ® ® ■ ■ ■ ® Ifc ® ■ • ' ® ^n- 

fe=l 

By [2], Proposition 2.2, p(T) is an essentially self-adjoint operator. 
Pick £ € V, T an operator on V and A G R. Assume that 

T is essentially self-adjoint, the vectors {T k ^}^L belong 
to its dense domain, span it, and are analytic for T. 

Define a t (0 = o(l [0 , t ) ® £), o*(0 = «*(l[o,t) ® 6, and p*(T) = p(l [0 , t) ® T). 
Then the corresponding q-Levy process is 



p t (t, T, A) = a$(0 + ot(e) + Pt (T) + At. 
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Let {X(t)} be such a process, let 2 be a subdivision of the interval [0,i), 
and let {Xi\ be the increments of this process corresponding to the subdi- 
vision intervals of I, X{ = X(di+\) — X(a{) for Ij = [a^, aj+i) E I. Then 

A k (t;X,I) = J2x? 

i 

and 

ii,i 2 ,...,i k 
distinct 

The stochastic measures A^(t;X) and ipk(t;X) are the limits of the above 
quantities as the size of the subdivision 6(1) tends to 0, if these limits exist. 
Similarly, for a /c-tuple of processes we can define 

A(t; (X W , X^ X^ ),1) = Y J X\ l) X\ 2) . . . Xi k) . 

i 

Such a fc-tuple form a multiple g-Levy process if they satisfy an extra com- 
patibility condition ([2], equation 1), similar to the one in (1). 

Lemma A.l ([2], Proposition 3.6). Let {X®(t) =Pt(&,2 1 i) A i )}£ =1 be a 
multiple q-Levy process. Then the q-cumulants 

R(t;(X^,X^ 2 \...,X^)) = lim E[A(t; (X^^X®, . . . ,1^),!)] 

<5(X)— >0 

are well defined, and equal to 

R(t; (X m ,X i2 \ . . .,X (t) )) = | t^ft, J] Tj&y ',f k >! 

Proposition A. 2. Let X(t) = pt(£,,T, \) be a general one- dimensional 
q-Levy process. Then the limit defining A^faX) exists in the L 2 -norm with 
respect to E[-], and equals 

Y(t)= Pt (T k - 1 t,T k ,(Z,T k - 2 0). 

Proof. Condition (1) implies that any A:-tuple of processes whose com- 
ponents are X and Y is also compatible. It suffices to show that 

Um o ((A fc (t;A' ) T) - Y(t)) 2 tt,tt) = 0. 
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First expand 



= E *f + E *W - E * - E ^4 - E xft - E ^ + Y " 

= A 2k (t; X,l)+J2 X k X k -A(t;(Y,X,...,X),l)-J2 
- A(t; (X, . . . , X, Y),I) - E ^ + y(t) 2 . 

From the pyramidal independence of the increments ([2], Lemma 3.3), it 
follows that 



(2) 



E 



E^W 

i+3 



E E[Xf ]E[X*] = R k (t; X,l) 2 - E E[Af ] 2 . 

i^j i 

By combining Proposition 5.2 with Lemma A.l and the moment-cumulant 
formula ([2], equation 3) 

E[X(t) k ] = E 

it follows that ELY(i) fc ] = J2neP(k) t^Rw^) = 0{t). Since the increments of 
the process are stationary, we also know that E[Xf] = ELY(|/j|) fc ]. Therefore, 

lim VE[X fc l 2 <C lim V llA 2 < C lim t6(I) = 0. 
6(1)^0^ <S(X)->0^ ~ 5(J)-0 

We conclude that the limit of the expression (2) is R k (t;X) 2 . Similarly, 

E[Y(t)]R k (t;X). 



lim E 

<S(X)->0 



E^4 



lim E 

<5(Z)->0 



Therefore, 



lim ((A fc (t;X,J)-y(t)rn,fi) 
5(X)->0 

= i? 2 fc(t; X) + R k {t; Xf - R(t; (Y,X, . . . ,X)) — E[Y(t)]R k (t; X) 

- R(t; (X,...,X, Y)) - R k (t; X)E[Y(t)} + E[Y(t) 2 ] 

= (£,T 2k ~ 2 + {i,T k ~ 2 i) 2 - {T k - 1 Z,T k - x £) - (t,T k - 2 £}(Z,T k - 2 £) 

- (^,T k - 1 T k - 1 ^) - (C,T k ~ 2 0(^T k ~ 2 
+ (T k - 1 ^,T k - 1 ^) + (C,T k ~ 2 2 

= 0. 
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□ 

PROOF of Proposition 5.3. For the g-Brownian motion, T = 0, while 
for the g-Poisson process, T = Id. So the result follows from the preceding 
proposition. □ 

Proof of Proposition 5.5. It suffices to show that the stochastic 
measures satisfy the same recursion relations as the orthogonal polynomials 
Cq >n . That is, we will show that 

lim \\X(t)Mt;X,l)-iP n+1 (t;X,l) 
(3) «(z)-*o 

-[n] q ij n (t;X,l)-t[n] q ^ n . 1 (t;X,l)\\ 2 = 0. 
Indeed, if that is the case, then 

Vwi (t; X, 1) = X(t)^ n {t; X,X) + [n] q i/> n (t; X, 1) + i[n]^ n _i(i; X, 1) + A(2) , 

with I? — lim^x)— »o = 0. By induction, the right-hand side converges 
in L 2 to 

X(t)rp n (t;X) + [n] q rp n (t;X)+t[n] q rp n - 1 (t;X) 

[for the first term, we use an argument similar to the proof of (1)]. So the 
left-hand side also converges in 1? . Moreover, also by induction, the limit of 
the right-hand side is 

X(t)C q , n (X(t)) + [n] q C q>n (X(t)) +t[n] q C q>n - 1 (X{t)) = C 9 , n+1 (X(i)). 

We will omit X,I and t in the notation. Expanding the norm in (3), we 
get 

E[(ip n+1 + [n] q ^ n + [n) q til) n - 1 - Xij) n ) 
x (ip n+1 + [n] q ip n + [n] q tip n -i - ip n X)\ 

= E[^„+l^n+l] + WgE^n^n] 



(4) 



+ [n]^ 2 E[^ n _i^ n _i] + E[X^ n X] 

- E[tp n+ iip n X] - E[Xtl) n il) n+ i] - [n} q E[ip n ip ri X} 

- [nlqElXipnipn] - [n] q tE [ip n _ i ip n X] - [n] q tE[Xip n iJj n ^ 1 ] . 



Combining the general linearization formula (2) and the specific form (6) of 
the cumulants of the q-Poisson process, 

It is easy to show by induction that 
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This follows from the fact that 1 is connected by tt to exactly one element 
2n — k, their class crosses exactly k — 1 other classes, and J2k=i Q k ~ 1 = Mg- 
We conclude that 

E[M n } = [n] q \t n . 

We similarly simplify the other expressions in the sum (4) . We treat in detail 
the most complicated term 

E[Xtf; n ^ n X] = iE[^„] + [ngE[^ n ] + (1 + g)[n]^ 2 E[Vn-i^„-i] 

= [n] q \t n+l + [n]\[n] q \t n + (1 + q)[n] q [n} q \t n+1 . 

Denote by V the join of set partitions. In the first three terms in the sum 
below, 7T is a partition in Viin, n) induced on the subset {2, . . . , In + 1} of 
{1, . . . ,2n + 2}. In the last two terms, tt is a partition in Viin ~ — 1) 
induced on the subset {2, . . . , 1 + k, . . . ,2n + 2 — j, . . . , 2n + l} of {1, . . . , 2n + 
2}. Using (1), 

= 5Z St 7rV {( li2n+ 2),(2),...,(2n+l)} 
7r[{2,..,2n+l} 

+ 5Z 

OT W{(l,l+fc),(2n+2-j,2n+2),(2),...,(2n+l)} 

■7TGV2 (n,n) k,j=l,...,n 
7rt{2,...,2n+l} k~j modTr 

+ 5Z H 

5T W{(l,l+fc),(2n+2-j,2n+2),(2),...,(2n+l)} 

■7TGV2 (n,n) k,j=l,...,n 
Tr\{2,...,2n+l}k^j modTr 

n 

+ 5Z 5Z St 

7rV{(l,l+fc),(2n+2-j,2n+2),(2),...,(2n+l)} 

k,j=l 7reP(n-l,n~l) 

7rt{2,...,r+fc,...,2n+2^j,...,2n+l} 

n 

+ 5Z 5Z St 7rV {( lj2n+ 2-j),(l + fc l 2n+2),(2),...,(2n+l)} 

fc,i=l TrePfn— l,n— 1) 

7rt{2,...,I+fc,...,2n+2^jr',...,2n+l} 

+ terms containing singletons. 

See Figure 3 for an illustration. Five types of partitions in it correspond to 
the five terms in the expression above. The classes containing the first (i.e., 
"1") and the last (i.e., "2n + 2") elements of the set are shown; the remaining 
classes consist of two elements each and are inhomogeneous with respect to 
the partition 7r n+ i in+ i. The dashed classes belong to the partition tt. Note 
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I (1+k) (2n+2-j) (2n+2) 

Fig. 3. Five types of inhomogeneous partitions without singletons obtained from the prod- 
uct V'lV'nV'nV'l- 



that in the third diagram, the dashed lines may cross each other as shown 
in the diagram, in which case the crossing is counted among the crossings of 
7T. Alternatively, the dashed lines may not cross, in which case one of them 
crosses a solid line, and this crossing is counted among the extra ones in 
the sum below. In the first diagram, there are no crossings between classes 
of 7r and the extra class. In the next three diagrams, the two extra classes 
cover (k — 1) and (j — 1) points, respectively. Because tt is inhomogeneous, 
each point covered by an extra class has to be connected to a point not 
covered by it, hence, extra {k — 1) + (j — 1) crossings are introduced. In the 
last diagram, the classes containing points {2, . . . , k} have to cross the extra 
class (1 + k, 2n + 2) , and the classes containing points {2n + 3 — j, . . . , In + 1} 
have to cross the extra class (1, 2n + 2 — j). In addition, the two extra classes 
also cross each other. So extra (k — 1) + (j — 1) + 1 crossings are introduced. 

Taking expectations, using the fact that the process is centered, and the 
specific form (6) of the cumulants, we obtain 



E 
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+ j2 ( E <z fe "¥ - y c(7r w e g fc -y"y cW t nN ) 

tt£V2 {n,n) V k,j=l,...,n k,j=l,...,n / 

mod7r fc^j mod7r 

n 

+ E E ^-y-y^tv- 1 

feJ=l7reP2(n-l,n-l) 
n 

+ E E iVV" 1 ^ -1 

feJ=l7rGP2(n-l,n-l) 

= [n] g !i" +1 + [ng[n],!t n + [n]> - l] q lt n+1 + g[ng[n - l],!*^ 1 . 
Similarly, 

E[i> n+1 il> n X] = E[X^ n+l ] = E^+iVn+i] = [n + l) q U n+1 , 
E[^ n X] = E[Xt/j n ^ n ] = [n] ff E[^ n ] = [n] q [n] q \t n 

and 

E[^ n ^ n X\ = E[X^ n -i] = [n],*E[^ n _i^ n _i] = 
Substituting these relations into (4), we obtain 

[n + lVr +1 + [n] 2 g [n] q H n + [n] 9 H ? !f + 1 

+ [n] q \t n+1 + [nf q [n] q \t n + (1 + g)N g [n] g !t n+1 
- 2[n + - 2[nf q [n] q \t n - 2[n] q [n} q \t n+1 

= [n] q H n+1 + q[n] q [n] q H n+1 - [n + 

= (l+g[n] g -[n+l] g )[n],!t" +1 =0. □ 

Proof of Proposition 2.2. Let F G A s - Then 

E[T^(t)]= hp E[Y^ n (i;X)]. 
5(1)— *0 

Since the limit exists, we may restrict Z to subdivisions containing s as an 
endpoint of one of the intervals. The above expression is a sum of terms of 
the form 

If I v u\ (£. [0, s] for some j, the corresponding X v ij\ is singleton independent 
from the rest of the terms in the product. Since the process is also centered, 
the resulting expectation is 0. As a result, 

E[Y^ n (t;l)]=E[Y^ ri {s;l)) 
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and so 

E[YMt)]=HYMs)}- 

Since this equality holds for an arbitrary Y E A s , we conclude that the 
conditional expectation of ip n (t) onto A s is ip n (s). □ 

A transition operator for a Markov process is called Feller if it maps Co (M) 
into itself. 

Corollary A. 3. Let C q ^ n be the scaled version of the continuous big 
q-Hermite polynomials and {X(t)} be the centered q-Poisson process: 

(a) Cq, n (X(t),t) is a martingale with respect to the filtration induced by 
the process {X(t)}, for every n. 

(b) Let 

oo 2 

H{X ' *' Z) = Jll + ztqK-(z q K/(l+zqK))(l-q)x- 

Then H(X(t),t,z) is a martingale. 

(c) The process {X(t)} is a Markov process with a Feller kernel. 

Proof. Let pt be the orthogonal projection from L 2 (]R + ,dx) onto the 
subspace L 2 ([0,t),dx). It can be extended to an operator on JT (? (L 2 (R_|_)). 
The conditional expectation onto A s is obtained by compression: 

Et[A] =ptApt. 

The increments of a (/-Levy process are pyramidally, and so singleton, inde- 
pendent. Thus, the first part of the corollary follows from Propositions 2.2 
and 5.5. It can also be obtained from the chaos decomposition property for 
the g-Poisson process, 

Cq t n{X(t),t)Q = 

The second part follows from the first one since 

oo -. 

H(x,t,z) = J2 TT -C q , n (x,t)z n . 
Note that the product defining H converges for all z since the sum 

converges. 
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The third part follows from the observations that the polynomials {C q ^ n } 
are, for every t, a basis for the polynomial ring, and polynomials are uni- 
formly dense in the space of continuous functions on the (compact) spectrum 
of X(t). Since the conditional expectation onto A s is norm-continuous, this 
implies that for any continuous /, it maps f(Xt) into the C*-algebra gen- 
erated by X(s). The existence of a Feller Markov kernel follows, see [5]. 
□ 
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